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Abstract 

In this paper, we propose some algorithms for the simulation of the distribution of certain 
diffusions conditioned on terminal point. We prove that the conditional distribution is absolutely 
continuous with respect to the distribution of another diffusion which is easy for simulation, and 
the formula for the density is given explicitly. 

1 Introduction 

The aim of this paper is to propose algorithms for the simulation of the distribution of a diffusion 

dxt = b(t, xt)dt + a(t, xt)dwt, xq = u, < t < T, 

conditioned on xt = v, where b and a are given functions with appropriate dimensions, and w is a 
standard Brownian motion. 

From the point of view of application, this allows to do posterior sampling when the diffusion 
is observed at instants {t\, ■ ■ ■ , t n } C [0, T]. 

Let us recall that in the usual conditioning (see, e.g., [7j), the distribution of the diffusion x 
conditioned on xt = v is the same as that of another diffusion y satisfying 

dy t = b(t,y t )dt + a(t,y t )dw t , y = u, < t < T, 

where 

b(t,x) = b(t,x) + [aa*]{t,x)V x {logp(t,x;T,v)), 

and p(s,u;t, z) is the density of x t ' u . However, this is not suitable for simulations because in 
general, one does not know the transition density p. 

We will prove that, in certain cases, the conditional distribution of the diffusion is absolutely 
continuous with respect to the distribution of another diffusion which is easy for simulation, and 
we give the explicit formula for the density. This leads to an efficient simulation algorithm. 

Two different cases will be considered: 



1. The matrix a(t,x) depends only on t, and b has the form b(t,x) = bo(t) + A(t)x + a(t)bi(t,x). 
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2. The matrix a(t, x) is uniformly invertible. 

This simulation algorithm can be applied to a problem of parameter estimation for a discretely 
observed diffusion as in jSJ. 

This paper is organized as follows: In Sect. 2, we recall a Girsanov theorem for unbounded drift 
which is essential for our simulation algorithm. In Sect. 3, we consider Case 1, and in Sect. 4, we 
consider Case 2. 

2 A Girsanov theorem for unbounded drifts 

This section is devoted to give a slightly generalized Girsanov theorem which will be used in the 
next section. We call a measurable function F(t, x) from M + x M. d to M. n locally Lipchitz with respect 
to x, if for any R > 0, there exists a constant Cr > 0, such that, for any (t, x, y) € K+ x W d x W 1 
with \x\ < R, \y\ < R, 

\F(t,x)-F(t,y)\<C R \x-y\. 

And on the metric space C([0, T]; W 71 ), we define the filtration {J~t}t to be the natural filtration of 
the coordinate process. 

Theorem 1 Letb(t,x), h(t,x), a(t,x) be measurable functions from W + xM. d toW d , W 71 , andM. dxm 
which are locally Lipschitz with respect to x; consider the following stochastic differential equations: 



dx t = b(t,x t )dt + o-(t,x t )dw t , (1) 
dyt = (b(t,y t ) +a(t,y t )h(t,yt))dt + a(t,y t )dwt, yo = x , (2) 

on the finite interval [0, T]. We assume the existence of strong solution for each equation. We 
assume in addition that h is bounded on compact sets. Then the Girsanov formula holds: for any 
non-negative Borel function f(x,w) defined on C([0, T]; x C([0, T]; W 71 ), one has 

E y [f(y,w h )} = ^.[/(x,u;)ei'o T/l *^) d?i "-5/o T IMM t )l 2 d i]) (3) 
E x [f{x,w)] = E y [f(y,w h )e- So >(t.ift)d««-§ fo ' \h(t,y t )\ 2 dt^ ^ 

where w\ = u>t + f h(s, y s )ds, and h* stands for the transpose of h. 

Proof: We assume first that the positive supermartingale 

M t = exp|y h*(s, x s )dw s - ^ J \h(s, x s )\ 2 ds 

is a martingale under P x which will be proved later. In this case wt = wt — Jq h(s, x s )ds is a 
Brownian motion under MtP x , leading to a solution (x,w) of (0): 

dx t = b(t, x t )dt + a(t, x t )h(t, x t )dt + a(t, x t )dw t . 

As b(t,x),h(t,x),a(t,x) are locally Lipschitz with respect to x, pathwise uniqueness holds for 
Q and ©• The standard Girsanov theorem implies that (j3J) holds. 
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We prove now that Mt is a martingale. For any R > 0, consider the stopping time 

tr = inf{t > : \x t \ >R} AT. 

Taking into consideration that h is locally bounded, we have, according to the Girsanov theorem 
for bounded drift: 

P y\T TR = M r R P x \r TR - 

Hence 

E X [M T ] > E X [1 TR=T M T ] = E X [1 TR=T M TR ] = P y [T R = T] 

which converges to 1 as R — > oo. It implies that E x [Mt] = 1, and M is a martingale. 

Finally, follows in the same way. I 

3 Case when a is independent of x 

We assume here that xt has the specific form 

dx t = (a t h(t,x t ) + A t x t + b t )dt + a t dw t , x Q = u, (5) 

where at and At are time dependent deterministic matrices and h(t,x), bt are vector valued with 
appropriate dimension. 

For example the 2-dimensional process (x, y) which satisfies the following SDE: 

dx t = y t dt (6) 
dy t = b(t, x t ,y t )dt + adw t (7) 

and which is the noisy version of it = b(t,xt,xt), see, e.g. pQ. 
We shall prove the following result: 

Theorem 2 Assume that At, bt and at are bounded measurable functions oft with values in M. dxd , 
M. d and M. dxm , respectively. Assume also that h(t,x) is locally Lipschitz with respect to x uniformly 
with respect to t with values in W 71 , and locally bounded; and the SDE^) has a strong solution. 
Moreover, we assume that a admits a measurable left inverse almost everywhere 1 , denoted by a + ; 
and that h,A,b and a + are left continuous with respect to t. Then, 

(i) the covariance matrix R s t of the Gaussian process £t corresponding to with h = is given 

by: 

Rst = Ps / P~ l a u a* u P-*du P*, 



Jo 

where 

dP 

— l - = A t Pt, Po = Id, 
dt 

andP-* = (p- 1 Y; 



lr This requires essentially that a* a is almost everywhere > 



3 



(ii) the distribution of the process 

Pt=tt~ RtTR$ T (tT - V) (8) 

is the same as the distribution of £ conditioned on £r = v (M + stands for the left pseudo-inverse 2 
of M). For any nonnegative measurable function f, 

E[f(x)\x = U,X T = V}= CE r/(p) e /o T ^(*.POK + ^-^ + (A t p t +fe t )di)-|/o T |IMt,POI| 2 dtj ^ (g) 

where C is a constant depending on u, v and T. 

Proof: (i) The formula for R st is classic and comes from ^ = P t / * P~ l (b u du + a u dw u ) + -Pt£o> 
see e.g. jSj. 

(ii) Let us first recall that if (Y, Z) is a Gaussian vector, the distribution of Y conditioned on 
Z = zq coincides with the distribution of another Gaussian vector Y — RyzR\ z {Z — z o)> where 
R% z is the left pseudo-inverse of Rzz\ its covariance is Ryy — RyzR zz RzY- Taking Y as the 
vector (£ tl , • • • , £t fc ), and Z = £t, we observe that, defining the process p by (JHJ), (p^ 1 , ■ ■ ■ ,Pt k ) has 
the same distribution as that of (^t 1 , • • • conditioned on £t = v . And the covariance of pt is 

Cst = -Rsi — RsTRxT^Tt- 

Denote by p\ the process in particular for any nonnegative measurable function </?(•), 
^■[v(C)] = / -^[yO^^MTCdv) where /iy is the distribution of £r. For any nonnegative measurable 
functions / and g, 

E[f(x)g(x T )} = J E[/(0 5 (e r )ei'o T?t *('^) du "-5/o T ll^)ll 2 d i] 

= E[f{i)g{i T )e^ ^^^t^t-^t^t+b^-l^ ||ft(t,6)|] 2 dt]_ ( 1Q ) 

Given a sequence of partitions (A n ) n >i of [0, T\: 

A n = {t%<q<---<q n =T} 

with |A n | = max <j<fc n _i(t™ +1 — i") — ► 0, and a continuous stochastic process X, we define: 

k n 1 

S n (X) = h*(t2,X t n)a+(X %i -X t n). 

i=0 

Then 

E[\s n (0 - s m (0\ ai]= / £[|s n (p<0-s m (^)| Ai] / u T (d v ), 

which implies that S n (p v ) converges in probability fix- Hence, we can define J Q T h* (t,p^)a^ dp^ 
as the limit (in probability P(g> fix) of the sequence S n (p v ). Obviously, this limit is independent of 
the sequence of partitions (A n ) n which satisfies |A n | — » 0. 

Finally, defining the continuous function Qn(x) = N A x, x > 0, we have 

1 rT 



E[e N (f(£)g(£ T )e Sn ®-fo ^(«0^ + (^6+^)dt-|/ J ||fc(t,£ t )|| 2 dt)] 



2 M + = (M* M) M* and the symmetric matrix is inverted by diagonalisation with 1/0 = 



4 



Taking the limit first in n and then in N, and returning to (|10[). we deduce: 

E[f{x)g(x T )} = [ E[f(p v )efo ^ ^W)( CT t +d P?- CT t + ( A ^+ fe Odi)-| J T 11^(^)11^^^)^ ( du ) 



which implies Q and C is the value of the density of \xt with respect to the distribution of xt at 
v. I 

As the Brownian bridge, we have: 

Proposition 3 Let us assume that Mt = §7 P^ 1 a u o'^ l Pu*du is positive definite for any t £ [0, T). 
T/ien the distribution of the process p is the same as that of q which is the solution to the following 
linear SDE 

dq t = A t q t dt + b t dt + <J t a* t P t -*Mt l (P^(E[t t ) - q t ) - P£\E[Z T ] - v))dt + a t dw t , (11) 
with qo = u. 

Proof: The matrix Q t = PtM t is solution to Q t = (A t - a t o* t Pf* Mf 1 Pf^Qt, implying that the 
covariance of qt can be rewritten as follows: for s < t, 

Q s I'* Q-^uolQ^duQ* = Q s f M^P^a^P^M^duQ: 
Jo J 

= Q s (M' 1 - Mq 1 ) Q* t 

= P S M S (M~ 1 — M Q ~ 1 )M t Pt 

= P s {M Q -M s ){Id-M^{M Q -M t ))Pt 

= C s t- 

On the other hand, from ©, the expectation pt of the process pt satisfies 

^.pt-Atpt-bt = -a t a;P t -*PTR T T( E M-v). 

Elementary algebra shows P^R^t = —Q^ 1 (RtTRj'T ~ -^t^T 1 )' hence 

^p t - A t pt - h = o- t o* t P t -*Qt\R tT R T l T - P t P^){E[i T ] - v) 
= a t a* t P-*Qf{E%] - p t - P t P~\E[i T ] - v)) 

which is the equation satisfied by E[qt\. The conclusion follows by noting that both p and q are 
Gaussian processes. I 

Remark. Mt is positive definite for any t £ [0, T) if and only if the pair of functions (^4, a) is 
controllable on [t, T] for any t £ [0,T). See, e.g. (> for some discussions. 

Example. Consider the 2-dimensional stochastic differential equation defined by ()6I7|) . where 
a 0. Let us assume that b is locally Lipschitz with respect to (x,y), and this equation admits a 
strong solution (the strong solution exists if there exists a Lyapunov function, see, e.g. pQ). Then 
we have: 

E[f(x,y)\(xo,vo) = u, (x T ,y T ) =v]=CE r/ ( p, g ) e -- 2 /o T ^ t , 9t )d^-^/ ( r^ t , qt ) 2 d*l (12) 
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where (p, q) is the following bridge starting from (po, go) = u: 



Pt \ _ ( zt \ _± ( t(3T-2t) -tT(T-t) \ ( zt- Vi \ 
Qt J \z t ) T*\ 6(T-t) T(3t-2T) J\z T -v 2 )> { ) 

with zt = u\ + tu2 + cr / w s ds, if = U2 + cru^; 

■/ 

or (p, q) can be chosen as: 
dp t = q t dt, 

^-^r)*^ <14> 
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4 a invertible, general b 



4.1 Bounded drift 

Let us consider the following SDEs: 

dx t = b(t,x t )dt + a(t,x t )dw t , x = u, (15) 
dyt = b(t,y t )dt- ^— - dt + a(t,yt)dw t , y = u. (16) 

Remark. If b = 0, and a = Id, then x is a Brownian motion. It is well known (see, e.g. |0j) 
that the law of the Brownian motion x conditioned on xt = v is the same as that of the Brownian 
bridge y satisfying the following SDE: 

dyt = ~ YZjdt + dwt, yo = u. 

The form of SDE (|16j) is inspired by the above SDE in order to fit the simplest case: the Brownian 
bridge case. 

The objective of this section is to prove that the distribution of x (solution of (|15[l) conditioned 
on xt = v is absolutely continuous with respect to y (solution of (1161) ) with an explicit density. We 
shall assume some regularity conditions on b and a here. 

Assumption 4.1 The functions b(t,x) and a(t,x) are C 1,2 with values in M. d and M. dxd respec- 
tively; and the functions b, a, together with their derivatives, are bounded. Moreover, a is invertible 
with a bounded inverse. 

Let x s ' u be the solution of ()15|) starting at s £ Under Assumption 4.1, x is a strong 

Markov process with positive transition density. For (s,u) G [0, T], we denote p(s,u;t,z) to be 
the density of x s t ,u . Then there exist constants m, X, M, A > 0, such that the density function 
p(s, u; t, z) satisfies Aronson's estimation [2] : for t > s, 

d A 1 2 — d A| z — 

m(t — s)~ze < p(s, u; t, z) < M(t — s)~^e *- s . 

We first study SDE (|TH|). 

Lemma 4 Let Assumption 4-1 hold. Then the SDE A16\) admits a unique solution on [0, T). More- 
over, \\m. t ^Tyt = v,a.s. and \yt — v\ 2 < C(T — t) log log [(T — t)~ x + e],a.s., where C is a positive 
random variable. 



Proof: The fact that the SDE (|16j) admits a unique solution on [0, T) is classic. Applying Ito's 
formula to ^Er > we deduce easily the following: 



yt — v u — v 
T - t ~ T 



+ f {T-s)- 1 b(s,y s )ds+ I (T-s)- l a{s,y s )dw s . 
Jo Jo 



For each i, f Q (T - s) 1 a(s,y s )dw s j ,t > 0} = {^2 j=1 J (T - s) 1 a ij (s,y s )dwi,t > 0} is a 
continuous local martingale, and its quadratic variation process Tt = J Ylj=i(T ~ s )~ 2<7 ij( s i Vs)ds 
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satisfies Tt — > oo as t — > T, and < for a constant c > 0. Applying Dambis-Dubins-Schwarz's 
theorem, for each i, there exists a standard one-dimensional Brownian motion B l , such that 



( jf (r-a)- 1 <r(«,y.)dw a ). = B\r t ),t> 0. 



Taking into consideration of the law of the iterated logarithm for the Brownian motion B l , the 
conclusion follows easily. I 
Now we can state the main theorem of this section. 

Theorem 5 Let Assumption J^.l hold. Then 

E[f(x)\x T = v] (17) 
rT 2y}A t (y t )b t (y t )dt + y* t {dA t {y t ))y t + £. . d(A?(y t ), y\y>) ' 



CE 



/(y)exp 



2(T - t) 



where A(t,y) = (a(t,y)*) l a{t,y) 1 , yt = yt~v, and (-, •) is the quadratic variation of semimartin- 
gales. 

Remark. Prom Lemma 4, the integral in ((17)) is well defined. 
Proof: Let fix) be an ^-measurable nonnegative function, t < T, then 



E[f(y)] = e 



/(x)exp 



1 {aj\x s )(x s -v)Y A 1 /■*.(77 1 ( a ; a )(x a -t;) 



T-s 



-dw 



T-s 



l ds 



(18) 



On the other hand, Ito's formula gives: 



d \\o- l {t,x t )(x t -v)\\ 2 = 2 (x t -v)*A{t,xt)dx t + \\a- l {t,xt){xt-v)\\ 2 At + d ■ dt 



T-t T-t {T-t) 2 

(x t -v)*(dA(t,xt))(x t -v) E 8J d(^(t,x t ),(^-^)(xj-^)) 
T-t T-t 
Combining the above equation with (|18j). we deduce that, 



T-t 



E[f(y)] 



CC t E 



fix) exp 



\a t l ix t )ix t -v) 
2(T - i) 



+ 



jx s - v)*A s (x s )b s jx s ) 
T-s 



ds 



T-s 



T-s 



where C > is a constant, and = (T — t) 2 . 
Or equivalently, 



E[f(y)<p t ] = CC t E 



f{x) exp 



\ajt,x t ) 1 jx t - v) 
2(T - t) 



where 



ft = exp 



* y* s A s (y s )b s (y s ) j l /■* y:(dvl s (y s ))y s , d(A^(y 8 ), y\yl 



T-s 



-ds 



T-s 



+ 



T-s 



(19) 



(20) 
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Note that {<ft,t £ [0, T]} is a well defined continuous process, thanks to Lemma|l] 
Putting / = 1 in ffTTTf) . we deduce then: 



E[/(y)^] 




/(x)exp j 


Hcr^iEi)- 1 ^-!;)!! 2 V 
2(T-t) JJ 


EN 




exp {- 


llo-fecct)- 1 ^*-!;)!! 2 !] 
2(T-t) / 





Assuming that /(x) takes the form f(x) = g(x tl , ■ ■ ■ , x tN ), < ti < t 2 < • • • 
<7 6 C&(R ), and letting i — > T, from the Lemmas [7| and |H] in the Appendix, we get: 



(21) 

< t N < T, 



E[/(y)yr] 

This completes the proof of the theorem. 



E[f(x)\x T = v}. 



Remark. For practical implementation, it is useful to note that the second and third terms of 
the integral in (jTZJ) are the limit of £ $ k (A(t k , yt k ) - A(t k - 1 ,y tk _ 1 ))y tk 2 (T-t k ) ■ 

4.2 Unbounded drift 

Let us now consider the following SDE: 

dx t = b{t, x t )dt + a(t, x t )dw t , xq = u, (22) 

where the drift b can be unbounded. We assume instead 

Assumption 4.2 The function a(t,x) is C 1,2 with values in R rfxd ; the function a together 
with its derivatives are bounded; and a is invertible with a bounded inverse. The function b is 
locally Lipschitz with respect to x and is locally bounded. Moreover, the SDE (|2*2*)l admits a strong 
solution. 

Combining the Theorems ^ and |51 we are able to prove the following 



Theorem 6 Let Assumption hold, and y be the solution of 

yt-v 



dy t 
Then, 



t 



dt + a(t, yt)dw t , y = u. 



(23) 



E[f(x)\x T = v] 



CE 



+ 



T y*(dA(t,y t ))y t + £\, d(A^(t, y t ), y\y> 



/(y)exp j 
T (b*A)(t,y t )dy t 



o 



2(T - t) 
\ J T \a- l b\ 2 {t,y t )dt 



where A(t,y) = a(t,y) *a(t,y) , yt = yt — v, and (-,-) is the quadratic variation of semimartin- 
gales. 
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Proof: Let x be the solution of: 

dx t = a(t,x t )dw t , xo = u. (24) 
Then, from Theorem ^ for nonnegative measurable functions / and g, 

E[f(x)g{x T )} = E[f{x)g{x T )e^^ A )^)^-\So |o- 1 6| a (*^*)dtj 



It remains to apply Theorem El I 
Remark. If the drift b is bounded, both formulas in Theorems [5] and are available. Unfortu- 
nately, it is difficult to compare the efficiency of simulation when applying these two formulas. 
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5 Appendix 

Lemma 7 Let < t\ < t 2 < ■ ■ ■ < t N < T, and g £ Cb(WL Nd ). Then, putting 



ipt = exp 



2{T -t) 



lim ~^^§tu ' XtN ^ = E ,xt 2 ,---, xt N )\xT = v]. 
PROOF: For any t E (ijv,T), 

E[g(x tl ,x t2 ,- ■ ■ ,x tN )ipt] 



E 

where 





i $ 9 (t,z)exp j 


f llaCM)" 1 ^)!] 2 








i $i(t,z)exp <j 


r n<7(t, Z )-i(^)ii 2 




L 2(T-i) 



* 9 



(t,z)= / g(zi, ■ ■ ■ ,z N )p(0,u;ti,zi) ■ ■ ■ p(t N , z N ;t, z)dzi ■ ■ -dz N 

jR Nd 

which is continuous thanks to Aronson's estimation. Evidently, z) = p(0,u;t, 
Moreover, applying a simple change of variable z = v + (T — t)^z', 



\a(t, z) l (z — v) 



1 2 



= I .,(*,. + ( t - t )^',exp {- H^±(r_t)M^T j dz , 

..(T,.) /, a exp{- »^f V » 2 }d,'. 

Hence, 

lim E[g(x tl ,x f2 ,---,a; fjv )^] = gg(7>) 
tir Efofe] *i(T,«)' 

from which we deduce (|25jl by the Bayes formula, since 

$ g (T,v) = / • • • ,z N )q(z 1} ■ ■ ■ ,z N ,v)dzi ■ ■ -dz N , 

where q is the density of (xt 1 , • • • , xt N ,xt)- 
Lemma 8 

lim E[\<p t - <pt\] = 0. 
We need the following two propositions to prove this lemma. 
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Proposition 9 (i) There exist two constants c\ > 0, C2 > 0, such that 

c 1 <C t E[^ t ]<c 2 ,yte[0,T), 

where 



C t = (T-t)-2. 
(ii) There exists a constant C3 > 0, such that 

E[(pt] <c 3 , Vi G [0,T). 

PROOF: (i) We note that 

C t E[ip t ] = (T- t)~* J exp|-Ml|l_^£_^| p( o,n;t,z)dz. 

We get easily the conclusion taking into consideration of Aronson's estimation after a change of 
variable z = v + (T — t)^z' '. 

(ii) It follows from (JT£J) and (i). I 

Proposition 10 For any e > 0, there exists an adapted bounded process at such that 

dC t A = dM t + atiCt^t) 1 - 6 ^ - t)~ h dt, h = £ -^- 

where (Mt)o<t<T is a, martingale. 

Proof: Set xt = xt — v, pt = \\a (t, x t )x t \\, and A t = o~~*{t, xt)o~~ l (t, xt). We have 

d=£- = 2^^ + 7 ^dt+^dt+^^ + ^Vd(^^ 
T-t T-t (T-t) 2 T-t T-t T _ t Z^\<' 



x* t a(t,x t ) *dw t p 2 t d Pt+Pt A . , Pt 

= 2 — h 7— rrrdt + — at + r t — at + — 

T-t {T-t) 2 T-t T-t T-t 

where r t and r' t are two adapted bounded processes. Hence we get: 



2(T-t) 2 r \T-tJ 8 1 \T-t 

si 

where r" is an adapted bounded process. For any e > 0, e 6 2 \x\ ,k = 1,2,3,4, are all bounded 
functions, then there exists a constant c £ > such that 

P 2 t+Pt . Pt+Pt 



T-t (T-t) 2 J ~ y/T^t' 
Hence, 

dCtiJt = dM t + (C t ilH) 1 - e (T-t)- h r%'c e dt, 
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where r' t " is still an adapted bounded process. 

Let us now return to the proof of Lemma |H1 
Proof: First, from Fatou's lemma and Proposition 

E[(Pt] < limmf -E^t] < C3. 
We choose to £ (0, T) which is close enough to T, and A large enough, and put 



1 , _ o . _,_ ... A 



a = inf{t < t < T, C t ip t < il = mf{to <t<T,p 2 t > 2(T - i) log ^^}. 

^ (T-t)5 



Under the distribution of x , <r < T a.s. However under the distribution of y,, limyi_> +00 a = T, 
a.s., taking into consideration of Lemma 01 We have, from (|21|) . 

E[ip t l u<t ] E[jj t l a<t ] 1 
E[(f t \ E[ip t \ ci 

On the other hand, from Proposition EH with a fixed e G (0, 1/d), 

dCtipt = dM t + atiCt^tf-^T - t)~ h dt, 

i.e., ^ 

C t ^ = + M t - M a + f a s {C s ^ s ) l - e (T-s)- h ds. 

J a 

Hence, 

E\C t ^ t \ a<t \ <A~ x + a\ E[C s ip s l a<s ] l - £ {T - s)~ h ds, with a = sup IKIU. 

Ji t 

Therefore, E[Ctipt^a<t] is bounded by ut which is the solution of the following differential 
equation, 

du t = au]~ e {T - t)~ h dt, ut = A" 1 ; 
and this equation has an explicit solution: 



[■^[(T ~ to) 1 -* 1 - (T - t) l ~ h ] +A-^ < {c (T - to) 1 -* 1 + A~ £ } ] 



where Co > is a constant. We get finally, 

E[tp t lt<cr] _ 1 E[cp t la<t] 



E[ft] E[pt] 

> l--( Co (T-t ) 1 ~ h + A- £ ) 1 / £ . 
ci 

We note that {(ptlt<a}t is a uniformly integrable family due to Novikov's lemma, since we have 

<r f f tAa \y s -y\ 2 1 1 r CT l^-^l 4 , ,2, \ 

lt<<r<Pt < C exp <J T- s ~ 2 J \T - s\ 2 J ' 
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where for fixed A, C is a positive constant and vt is an adapted bounded process. 
Taking the liminft^T, we get, 

v E[LPTl -r] . > l-i(co(T-t )^ + A-^. 
lim sup t ^ T E [ipt\ c\ 

Since l a =T converges to one a.s. as A — > oo, we get 
^ > l--(c (T-to) 1 -' l ) 1 / £ . 



lim sup t _> r £7[vt] ci 
It remains to let to — ► T to get: 

lim sup E[(p t ] < E[(px]- 
t-*T 

Hence, 

lim E[p t ] = E[tp T ], 
and we finish the proof by Scheffe's lemma (see, e.g. 4 j). 
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